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ABSTRACT 

Using a semi-analytical approach we investigate the characteristics of predictions for 
the masses and metallicities of the baryonic matter in and around galaxies made by 
three galaxy formation models. These models represent three different feedback sce¬ 
narios: one model with purely ejective feedback, one model with ejective feedback 
with reincorporation of ejected gas, and one preventative model. We find that, when 
the model parameters are adjusted to predict the correct stellar masses for a range 
of halo masses between 10^° to 10^^ M 0 , these three scenarios have very different 
predictions for the masses and metallicities of the interstellar and circum-galactic me¬ 
dia. Compared with current observational data, the model implementing preventative 
feedback has a large freedom to match a broad range of observational data, while the 
ejective models have difficulties to match a number of observational constraints simul¬ 
taneously, independent of how the ejection and reincorporation are implemented. Our 
results suggest that the feedback process which regulates the amounts of stars and 
cold gas in low-mass galaxies is preventative in nature. 

Key words: galaxies: evolution; galaxies: formation; galaxies:ISM; galaxies: stellar 
content 


1 INTRODUCTION 

One of the most fundamental questions in galaxy formation 
is how baryons are accreted into and expelled from dark 
matter halos as galaxies form. Using statistical approaches, 
a number of studies have revealed a picture in which the 
baryonic mass locked in to stars and cold gas is a strong 
function of galaxy mass (|Yang,_^fo_fc_jaiIlign_Bosch|l200j; 


Behroozi, Conroy fc Wechslej|2010r Moste r. Naab fc Whitel 

2OI3I : iPaoastergis et al.ll2012l : ILu et al.ll2014bl ). The Milky 


Way sized galaxies at the present day have the largest 
baryon mass fraction in stars, and this fraction decreases 
rapidly as the galaxy mass increases and decreases. How such 
a relation is established is one of the most important ques¬ 
tions in galaxy formation. A widely adopted galaxy forma¬ 
tion scenario attributes the low baryon mass fraction in low- 
mass halos to strong galacti c outflows generated by the feed¬ 
back of star formation (e .g. IPekel fc Silk|[l9^ : I Bower et al.l 
l200fil : ICroton et al.l[200^ . On the other hand, a number of 
authors have realized that preventing baryons from ever col- 
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I2OI5I) . Clearly, how feedback works in reality becomes a cen¬ 
tral question in galaxy formation theory. In this paper, we 
present a comparative study of three representative feedback 
models to understand the impact of ejective versus preven¬ 
tive feedback scenarios on their predictions for the masses 
and metallicities in different components of disk galaxies. 

The details of feedback processes are still poorly un¬ 
derstood. In the commonly adopted ejective scenario, feed¬ 
back of star formation is assumed to produce strong out¬ 
flows which can expel a large amount of gas mass from the 
host halo (e.g. iBenson et ahl l2003l : ISomerville et al.l I2OO8I : 
IGuo et aDl201ll ). However, the properties of the outflows, 
such as their masses and velocities, and how they are cou¬ 
pled to the halo gas are still uncertain. More recently, a 
more sophisticated scenario, in which ejected gas is allowed 
to reincorporate back t o the host halo after a certain time 
scale is proposed (e.g. iHenriaues et al.|l2013l ). In addition 
to the same uncertainties that the simple ejective model 
has, the ejective model with reincorporation also has un¬ 
certainties in the time scale for the gas reincorporation. The 
physics of these assumed processes is still poorly understood. 
In the preventive scenario, some early feedback is assumed to 
change the thermal state of the intergalactic medium around 
dark matter halos so that a fraction of baryons is prevented 
from collapsing into d ark matter halos in the firs t place. For 
example, as shown in lLu. Mo fc WechsleJ (l2015l) . when the 
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circum-halo medium is preheated by early feedback to a cer¬ 
tain level of entropy, the thermal pressure of the preheated 
gas can support the gas against gravitational pull, result¬ 
ing in a low bar yon mass fraction for low- mass halos. As 
demonstrated in lLu. Mo &: Wechslerl ll2015l ). such a model 
seems to be able to reproduce the observed stellar and cold 
gas mass fractions, sizes of disk galaxies, and star forma¬ 
tion histories of galaxies with mass comparable to or lower 
than the Milky Way, without invoking any strong outflow 
as adopted in the ejective scenario. An important question 
is then whether these different feedback scenarios can be 
distinguished with current observational data. 

Metallicities in different phases of a galaxy is expected 
to put interesting constraints on feedback models. Metals 
are produced by star formation and their distribution is 
affected by feedback. If, for example, the feedback drives 
weak outflows, metals would be found only in regions where 
star formation has taken place. On the other hand, if strong 
outflows are driven to large distances, metals produced by 
star formation can be redistributed not only in ISM (cold 
gas) but also in the hot halo gas and even in the intergalac- 
tic medium. In this paper, we a dopt a semi-analytic galaxy 
formation model introduced in ILu. Mo fc Wechsl^ (120151 ') 
to better understand how the masses and metallicities of 
baryonic matter in different phases are affected by different 
implementations of feedback. The goal is to examine how 
current and future observations may be used to distinguish 
between different feedback scenarios. 

The paper is organized as follows. In §(2] we describe 
the models adopted in this paper. We present our model 
predictions and comparisons with observations for the evo¬ 
lution of the baryonic masses in different components in §[3l 
and for the metallicities in different phases in §21 Our con¬ 
clusions and the discussion of the implications of the results 
are presented in §[5] Throughout the paper, we use a ACDM 
cosmology with Hm.o = 0.27, Ha.o = 0.73, Hb.o = 0.044, 
h = 0.70, n — 0.95, and ag = 0.82. 


2 MODELS 


We adopt the semi-analy tical model (SAM) introduced in 
ILu. Mo fc Wechslerl (l2015f l and focus on different implemen¬ 
tations of star formation feedback. To do this, we fix every 
other parts of the model, but vary the recipe for star for¬ 
mation feedback to study the impact of feedback on the ob¬ 
servational consequences of the masses and metallicities in 
stars, ISM and circum-galactic medium (CGM). We concen¬ 
trate on galaxies with mass comparable to or smaller than 
that of the Milky Way. We consider three representative 
models, which differ from each other by their assumptions 
on how feedback governs gas inflow and outflow. 

The f irst model is a typica l elect ive feedback model pre¬ 
sented in Mo fc Wechslerl (l2015l i. This model, referred 
to as Model-EJ in the following, captures the essence of 
many ejective feedback models in the literature in terms 
of how outflows driven by star formation fee dback affect 
the baryon contents of low-mass galaxies (e.g. Bovgr et al.l 


I2OO6I : iKang. Jing fc Sill] I2OO6I : ISomerville et al 


Ha¬ 
los are assumed to accrete baryons at a rate equal to the 
halo mass accretion rate multiplied by the universal baryon 
fraction fh = Hb,o/Hm,o, and a certain fraction of cold gas 


is assumed to be ejected out of the host halo as stars form. 
The ejected mass is assumed to be proportional to the star 
formation rate with a mass-loading factor that varies with 
the halo virial velocity I4ir. Based on the energy conser¬ 
vation argument, the mass-loading factor is assumed to be 
inversely proportional to 14 ^^ a-nd the mass ejection rate is 
written as 

A,. * 0.0 (52^) ", (1) 

where tp is the star formation rate (SFR), and old is the 
normalization of the mass-loading factor at K,ir = 200 km/s. 
The ejected mass is assumed to leave the halo and is added 
into an ‘ejected mass’ component. For the model considered 
here, we assume that the ejected mass never comes back 
to the halo to explore the maximum effect of the ejective 
feedback scenario. Based on these assumptions, the mass of 
the hot halo gas is given by 

Afhot — /bAfvir Afcold A/ej, (2) 


where Mvir is the virial mass of the halo, M* is the stel¬ 
lar mass and M^oid is the cold gas mass of the galaxy. 
The normalization factor, old, is tuned so that the 
model prediction matches the stellar mass of Milky Way 
sized halos as derived from recent empirical models that 
match galaxy stellar mass function s at multiple red shifts 
dBehroozi. Wechsler fc ConrovI l2013l : ILu et alTl2014al l. We 
find that old = 1 provides a reasonable fit to both the 
stellar mass and cold gas mass in Milky Way sized halos. 
Since the mass loading factor scales as lower-mass 

halos have larger mass-loading factors. Moreover, in order 
to explain the metal mass In the circum-galactic medium 
(CGM), we deposit a fraction, az, of the metal yield from 
star formation directly into the hot halo gas, and the rest 
is mixed uniformly into the ISM before being ejected out of 
the halo. We explore the effect of changing the fraction of 
this direct metal injection into the halo gas and found that 
az = 0.1 produces a good match to the OVI gas mass ob¬ 
served In C OS-Halos, assuming the OVI ionization fraction 
/oiv = 0.2 (iTumlinson et af] 1201 il l. A lower level of direct 
metal injection would predict a too low oxygen mass to be 
consistent with the observational constraints. 

The second model considered here is an extended ejec¬ 
tive model, in which the ejected gas mass is allowed to rein¬ 
corporate into the halo hot gas after a perio d of time. Specif¬ 
ically, we implement the model proposed bv iHenriaues et all 
(l2013ll . where the gas reincorporation rate is written as 


Mei = 


Mej 

treinc 


( 3 ) 


where treinc is th e time scale for the e jected gas to fall 
back into the halo. lHenriaues et al.l ll2013ll utilized a MGMG 
method to optimize a flexible version of the L-Galaxies SAM 
to match the evolution of galaxy luminosity function, and 
found the following model for treinc: 


treinc — O’ 


10“ Mq 

Mvir ’ 


( 4 ) 


where 7 is a free parameter tuned to be 7 = 1.8 x 10 ^'’ yr. 
To be consistent with the Henriques model, the prescription 
and parameter values of t he outflow model we adopt fol¬ 
low iHenriaues et al.l (l2013lf exactly. Specifically, a fraction 
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of cold gas is assumed to be heated and deposited into the 
hot halo gas reservoir and another fraction is assumed to be 
ejected and leave the halo. The ejected gas is then subject to 
reincorporation. The details of the model parameterization 
can be found in IGuo et al.l (1201 il l, and the detailed mod¬ 
eling of the reheating, ejection and reincorporation can be 
found in iHenriaues et ^ (l2013l ). In what follows, this model 
is referred to as Model-RI. Because this model already has 
reheated ISM to be mixed with the halo gas, we do not as¬ 
sume any additional channel of metal injection into the hot 
halo. 

The thi rd model is the preventativ e feedback model in¬ 
troduced in ILu. Mo fc Wechslerl (120151 '). In this model, the 
effect of feedback is not to eject gas out of halos but to pre¬ 
vent a fraction of the baryonic matter from collapsing into 
halos. This is achieved by assuming that the intergalactic 
medium (IGM) is pre-heated to a certain entropy level so 
that a fraction of baryonic mass will not be able to collapse 
into the shallow gravitational potential wells of low-mass ha¬ 
los owing to the thermal pressure of the preheat ed gas. As 
demon strated with hydrodynamical simulations in lLu fc Mol 
ll2007tl . a uniform entropy distribution of the IGM results in 
the fraction of baryon mass that collapses into a halo to be 
proportional to the halo mass squared, Mb oc M^ir, which 
is needed to match the shallow low-mass e nd slopes of loca l 
galaxy stellar mass and Hl-mass functions (IMo et al.ll2005ll . 
In addition, the halo gas in this model is assumed to estab¬ 
lish a h ot corona with an ex t ended central core in its density 
profile. ILu. Mo fc Wechslerl (l2015l 'l showed that this preven¬ 
tative scenario can reproduce the observed cold gas mass, 
stellar mass, star formation history, and size for galaxies 
with masses comparable to or smaller than that of the Milky 
Way without invoking any outflow after preheating. On the 
other hand, however, such a universal prevention cannot 
suppress star formation sufficiently for more massive galax¬ 
ies at late time. Thus, a “late quenching” process is needed 
for these galaxies. Here we consider a mechanism in which 
the stars and/or the central black hole in the central galaxy 
produce certain amounts of heating to off set the coo l ing o f 
the halo gas. Motivated by the results of |l^ fc Mol (120141') 
that the classical bulge mass and central black-hole mass in 
a galaxy with total stellar mass 10 ^°Mq < M* < 10^^M© 
both scale roughly as M*, we assume that the heating rate 
depends on M* as 

where the heating efficiency is set to be e = 5. In every time 
step, we calculate t h e halo cooling rate, C, as described in 
ILu. Mo fc Wechslerl ll2015h . and the heating rate. The net 
cooling rate is then 

Mcooi = maxp - H, 0] . (6) 

Moreover, we also allow a moderate level of outflow, assum¬ 
ing a constant mass-loading factor of unity for all galax¬ 
ies. This assump tion is consistent with directly observation s 
of outflows fe.g. iBouche et ^l2012l : iKacprzak et al.ll2014h . 
and is also allowed by statistical inference from an empirical 
model of galaxy evolution constrained by stell ar mass func¬ 
tions at different redshifts Mo fc Li3l2014^ . This model 
is referred to as Model-PR in the following. 

For all the three models, we follow not only the mass 


of bar yonic matter in each phase as in ILu. Mo fc Wechsi^ 
ll2015lj . but also metallicities in the different phases. For 
all the phases considered here, namely the halo gas, cold 
gas, ejected gas, and stars, we only track the average 
metallicity for each phase, but do not follow the details 
of the metallicity distribution. Throughout the paper, we 
adopt the solar metallicity as [12 - |- log(0/H)1© = 8.69 
(lAllende Prieto. Lambert&^^Agp hmd 2001: lAsplund et al.l 


l2009lj and Zrr, = 0.0134 (I Asplund et al.ll2009l) . and the total 
metal yield is taken to be p = 0.03, which is defined as the 
ratio between the mass of newly produced metals that are 
ejected into the ISM and the gas mass turned into stars. 
With the assumption that all galaxies have the same abun¬ 
dance ratio as the Sun, the oxygen yield is Poxy = 0.44p. 


3 BARYONIC MASSES IN DIFFERENT 
COMPONENTS 


To understand how the baryonic masses in stars, cold gas, 
hot halo gas and ejected components evolve with time in 
different models, we apply these models to typical halo mass 
accretion histories (MAHs) of two final halo masses, 10^^ M© 
and 10^^ M©. For simplic i ty, we adopt the fitting model of 
iMcBride. Fakhouri fc Mai (120091 '). 

Mhalo(^) = Mhalo.O (1 + 2)“ exp(-p2), (7) 


where the normalization Mbaio.o is the halo mass at 2 = 0 
and the parameters a and p determine the shape of the 
MAH. This functional form can be considered as a gener- 
alizati on of two other widely ad opted models proposed ear¬ 
lier bv lWechsler et al] (l2002|j and ivan den Bosm (l2002l l. and 
has been shown to accurately describe the assembly histories 
over a wide range of halo masses. It is also similar to other 
models proposed i n the literature (lZhao_eLalJ |200^|3, l2009l : 
IPekel et al.ll200^ . iMcBride. Fakhouri fc Ma ( 2009li showed 
that the parameters a and p have broad distributions. For 
our model, we choose the typical values of these two pa¬ 
rameters: a = 0.6 and p = 0.9 for Mhaio.o = lO'^^ M© 
halos, and a — 0.8 and p = 0.9 for = 10^^ M©. 

This model matches well the simulated MAH, as shown in 
ILu. Mo fc WechsI^ (l2015h . 

Following the spirit of SAM, we compute, at each time 
step, the halo gas accretion, hot gas distribution, radiative 
cooling, star formation and feedback, and keep track of the 
evolutions of the masses in stars, cold gas in the disk, hot 
gas in the halo, and the gas ejected out of the halo. The 
model predictions are shown in Figure[T] with the upper 
row for Mhaio.o = 10^^ M© halos, and the lower row for 
Mhaio.o = 10^^ M©. While the three models predict simi¬ 
lar albeit distinguishable histories for the stellar mass and 
cold gas mass, they predict very different histories for the 
halo gas and ejected gas masses. Model-EJ predicts the 
ejection of a large amount of gas and the ejected mass is 
kept out of the halo. Gonsequently the majority (more than 
half) of the baryonic mass originally associated with the 
halo is predicted to be ejected for both the halo masses. 
For Mhaio.o = 10^^ M©, the cold gas is the second dom¬ 
inating component, but for Mhaio.o = 10*^^ M©, hot halo 
gas mass is comparable to the mass of the stellar compo¬ 
nent while the cold gas mass is the smallest of all. The 
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Figure 1. The mass of different baryonic component as a function of time (redshift) for halos with a final mass of 10^^ Mq and 10^^ Mq at 
the present day. The predictions of different models are labeled by their names in different panels. The red and blue lines show the stellar 
mass and the cold gas mass in model galaxies, respectively. The orange line is the gas mass in the hot halo, and the purple line is the mass 
ejected out of the halo. The black dashed line plots the halo mass multiplied by the cosmic baryon mass fraction /h. The vertical bars 
mark the r anges of observat ional results. For stellar mass, the observational data c onstraint is taken f r om | Behroozi. Wechsler &: ConrovI 
ll2013lj and lLu et al.l For the cold gas, the observational constraint is fromlPapaa tergis ^ alj ||2012|~I. The c onstraint for the hot 

halo gas is adopted by bracketing data constraints obtained from iFang et al] 1 I 2 OO 6 I ') and iBregmaiTfc Lloyd Davies! ll2007l i. 


ejected mass never decreases with time because of the ab¬ 
sence of re-incorporation. The situation is very different in 
Model-RI. Although this model also implements a strong 
outflow, the ejected mass is re-incorporated back into the 
host halo at late times. Owing to the increasing reincor¬ 
poration timescale with decreasing halo mass, the ejected 
gas outside the halo is always the dominating fraction for 
low-mass halos, but it is a signihcantly lower fraction for 
A/haio.o = 10^^ M 0 . Furthermore, for the 10^^ Mq halos in 
this model, the ejected mass dominates the baryon mass 
budget only at very early time {z > 2), and this fraction 
decreases rapidly at low redshifts as the halo mass increases 
to the regime where re-incorporation becomes important. 
Finally, for Model-PR, the total baryon mass budget is gov¬ 
erned by the preheating entropy, which causes low-mass ha¬ 
los to be relatively poor in baryons. Since a low mass-loading 
factor of outflow is assumed in this model, the ejected mass 


simply follows the st ellar mass history with a delay caused 
stellar evolution (see ILu. Mo fc Wechslejl20lR for details). 
For Mhaio.o = 10^^ Mq, the cold gas mass is about 10® Mq 
at the present day, and the stellar mass is about half of the 
cold gas mass. For Mhaio.o = 10^® Mq, on the other hand, 
the stellar mass reaches about 3 x 10^® Mq while the cold 
gas mass is about one order of magnitude lower. For both 
Afhaio.o = 10^^ Mq and Mq, the hot halo gas mass is 
comparable to the stellar mass at the present time. 


We compare the predicted stellar and cold 
gas masses at the present day wit h the r es ults o f 
iBehroozi, Wechsler fc ConrovI (l2013lj , ILu et all (l2014ah 
and Papastergis et al. ( 2012ll . shown as the vertical red and 
blue bars. The observat io nal results for the hot halo gas ar e 
from iFa ng et all jl2006lj. iBregman fc Lloyd Davie j 120071), 


IWroet aLl l 2008f) . iTumlinson et al.l ll2011 1 and Gupt a et alJ 


1201 ^1 . and shown as the vertical yellow bars. Because 
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Figure 2. The rates of baryonic mass conversion among different phases as functions of time (redshift) for halos with a final mass of 
10^^ Mq and 10^^ Mq at the present day. The predictions of different models are labeled by their names in different panels. The red line 
is the star formation rate; the blue line is the radiative cooling rate of the hot halo; the purple line is the outflow rate. The black dashed 
line denotes the halo mass accretion rate multiplied b y th e universa l baryon fraction, q = 0.17. For comparison, the star formation 
rates inferred by iBehroozi, Wechsler fc ConrovI ll2013ll and iLu et aP ll2014a^ using empirical models are over-plotted as the green and 
orange bands, respectively. 


all the models are tuned to match the present-day stellar 
mass and cold gas mass of Milky-Way sized halos, it is 
not surprising that all models reproduce well the observed 
stellar and cold gas masses for Mq halos at z = 0. The 
measurements of hot halo gas mass have large uncertainties. 
Both Model-EJ and Model-PR can accommodate the 
observational results comfortably, while the prediction of 
Model-RI seems to be too high. 

For Mhaio.o = 10^^ Mq, both Model-EJ and Model-PR 
can achieve an excellent agreement with data for the final 
stellar mass, but Model-RI under-predicts the final stellar 
mass by a factor of 2 to 3. Even assuming strong outflows, 
both Model-EJ and Model-RI over-predict the cold gas mass 
by a factor of few. This over-prediction stems from the fact 
that a large amount of cold gas is settled in an extended 
gaseous disk where star formation is low due to the low 
surface density and feedback is not efficient to get rid of 
the cold gas. If we allow the gas disks to be smaller so that 


star formation can occur in a larger fraction of the disk, 
th e predicted gas mass can b e reduced. However, as shown 
in ILu. Mo fc Wechslerj (l2015li . to match the cold gas mass 
requires disk sizes to be much smaller than real disks. We 
also note that the overproduction of cold gas mass causes an 
underprediction of gas-phase metallicity (see next section for 
details). When the gas-phase metallicity is sufficiently low, 
the star formation rate surface density strongly depends on 
the gas surface density. This nonlinear dependence causes an 
upturn in stellar mass in the final 2 Gyrs for 10^^ Mq halo. 
In contrast to the ejective models, Model-PR can reproduce 
the stellar mass and cold gas mass for both 10*^^ Mq and 
10^^ Mq halos. In this model, the amount of cold gas settled 
in the outer parts of disks is smaller, as a large fraction of 
the gas is not able to make to the disk because of preheating. 

Figure [5] shows the rates at which the baryonic masses 
in different components transfer between each other. We 
focus on the coofing rate of halo gas, star formation rate 
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Figure 3. Cold gas mass - stellar mass ratio as a function of galaxy stellar mass at z = 0 predicted by d i fferent models, as labeled b y 
thei r names in different pane ls. The mean relations derived from observational data bv IPeeples &: Shankai] ll201lh . [Peeples et alj ll2014h . 
and iPapastergis et al.l ll2012li are shown by different lines as noted in the figure. 


of the cold gas, and the outflow rate as functions of time. 
For the outflow rate, we plot the rate at which the gas is 
leaving the galaxy. For Model-EJ, the outflow gas goes out 
of the halo, but for Model-RI a fraction of the outflow is 
retained in the halo. For simplicity, we do not distinguish 
between outflows going into the IGM and those going into 
the CGM. Since mass ejection is powered by star formation 
feedback in all the three models, the outflow rate is propor¬ 
tional to the SFR. The difference between models is in the 
assumed mass-loading factor. In Model-EJ, the mass-loading 
factor is proportional to and so low-mass halos have 

much higher outflow rate per star formation rate. Model- 
RI adopts a very large mass-loading factor to balance the 
gas re-incorporated into the gaseous halo, and so both the 
cooling rate and ejection rate are at least one order of mag¬ 
nitude higher than the SFR. In Model-PR, the cooling rate 
deviates significantly from the halo accretion rate because of 
preheating, and the star formation is largely controlled by 
the reduced rate of gas accretion. For comparison, we include 
the SFR as a function of ti me for both Mhaio.o = 10^^ Mq 
and 1 0^^ Mf7|, as inferred b y iBehroozi. Wechsler fc GonrovI 
||2013D an d inr et al.l (l2014al i from observational data. All the 
three models reproduce the average star formation history 
of 10^^ Mq halos quite well. For Mhaio.o = 10^^ Mq, Model- 
PR matches the empirical results well; Model-RI moderately 
under-predicts the SFR at low redshift; Model-EJ fails com¬ 
pletely: it predicts too high a SFR at high z and too low a 
SFR at low z. 

Figure [3] shows the present-day cold gas to stellar mass 
ratio /cold as a function of galaxy stellar mass predicted 
by the three models. The predictions are made by applying 
the models to a random set of halo mass accretion histo¬ 
ries s elected from the Bolshoi cosmological A^- body simula¬ 
tion llKlvpin. Trujillo-Gomez fc Primackll201ll l for the mass 
range between 10^° and 10^^ Mq. As one can see, the gas 
fractions predicted by all the three models for high mass 
halos are consistent with observational results. This again 
is because the model parameters are all tuned to match the 
properties of Milky Way sized galaxies. Both of the ejective 


models (Model-EJ and Model-RI) predict a steep decline of 
the gas fraction with stellar mass, so that both models sig¬ 
nificantly over-predict the gas fractions in low-mass galaxies. 
Varying the normalization of the mass-loading factor in the 
ejective models only shifts the amplitude of the predicted 
relation but has negligible effect on the slope of the rela¬ 
tion, and so cannot fix the mismatch. In contrast, Model- 
PR reproduces the observed /coia-M* relation quite well. 
This model predicts a large scatter in the cold gas fraction 
for high-mass galaxies. This is largely caused by the ‘late 
quenching’ we introduced, in which the suppression of gas 
cooling depends strongly on stellar mass. Also, star forma¬ 
tion in halos less massive than 10^^ Mq depends sensitively 
on halo mass accretion history in this model. Depending on 
whether a halo has acquired most of its mass before or after 
the onset of preheating, the total amount of stars that can 
form can be very different. This has the effect of increasing 
the scatter in the predicted fcoid-M, relation. 


4 METAL MASSES AND METALLICITIES IN 
DIFFERENT GAS COMPONENTS 

We also follow the histories of the metal mass in stars, cold 
gas, hot halo gas, and ejected gas of model galaxies hosted 
by 10^^ Mq and 10*^^ Mq halos, we adopt the usual instanta¬ 
neous recycling approximation for metals produced by start 
formation to be returned to ISM and CGM. The time delays 
between the nuclear synthesis of different heavy elements 
are ignored for s implicity. This appr oximation is supported 
by the results of iPeeples et al.l (l2014l l. who showed that the 
Oxygen mass, which is mainly produced by Type II SNe, is 
roughly a constant fraction of the total metal mass over the 
stellar mass range of our study. Figure 2] shows the histories 
of metal mass in each component as a function of time (red- 
shift) predicted by the three models. In each panel, we also 
show the total available metal mass, pM*/(l — R) (with p 
being the metal yield and R the return fraction), as star for¬ 
mation proceeds. We adopt R = 0.43, which is proper for a 
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Figure 4. The metal masses in different phases of baryonic matter as predicted by the three models (EJ, RI and PR, as indicated in 
the panels) for halos with final masses of Mq (lower panels) and 10^^ Mq (upper panels) at the present day. The red lines are tor 
the metal masses locked in stars, the blue lines for metals in cold gas, the orange lines for metals in gaseous halo, and the purple lines 
for metals that are ejected out of halo. The black dashed lines denote the metal mass available for the galaxy, which is approximated by 


Chabreir IMF. The ejective feedback in Model-EJ expels the 
enriched gas from the galaxy, and so most of the metal mass 
ends up in the ejected gas. In Model-RI, the ejected mass re¬ 
turns back to the halo after a short time delay for the Milky 
Way sized halo, but can stay in the halo for a longer period 
of time. Consequently, most of the metal mass in such a halo 
is in the hot halo gas. For 10^^ Mq halos, on the other hand, 
the ejected component contains a larger fraction of metals 
than the halo gas, because gas re-incorporation is slow for 
low-mass halos and much of the ejected gas stays outside of 
the halo. Furthermore, since low-mass halos are predicted to 
have a high cold gas fraction, the cold gas in such a halo is 
predicted to contain a large fraction of the total metal mass 
in Model-RI. In contrast, Model-PR predicts that the metal 
masses locked in stars and cold gas are always a substantial 
fraction of the total, while the metal mass in the halo gas is 
lower by a factor of a few. 

In Figure[5]we show the metal masses in different phases 
versus stellar mass for galaxies a.t z = 0 predicted by the 


three models, and compare the mode l predictions with th e 
compilation of observational data by IPeeples et al.l dZOldh . 
In addition, we also include a horizontal bar to represent 
the observational estimate of metal mass in hot halos de¬ 
rived from OVI absorptions of COS-Halos samples, assum¬ 
ing = 0.4 4p and the same ionizati on fraction /ovi = 0.2 
as adopted in iTumlinson et all llZOllh . All the models pre¬ 
dict that the metal mass increases with stellar mass, except 
for the ejected mass predicted by Model-RI. The general in¬ 
creasing trend is set by the total metal mass available, which 
is plotted as the black dashed line in the figure. Galaxies 
with higher stellar masses can produce more metals, but 
where the metals are located depends strongly on the feed¬ 
back model. 

In Model-EJ, a substantial fraction of the ISM is blown 
out of galaxies, and because all the ejected mass is assumed 
to leave the halo forever, most of the metal mass is predicted 
in the ejected component for galaxies of different masses. 
The model also has a small fraction of metals in the hot gas. 
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which stems from our assumption that a fraction of metal 
yield is deposited in the halo. For the particular parameter 
(az = 0 . 1 ) assumed here, the predicted metal mass in the 
hot halo is abont 2 x 10 ^ Mq for galaxies with final stellar 
masses ~ 10 ^ ° Mf?), which is consist ent with the observa¬ 
tional results llTumlinson et aLll201lli . The predicted metal 
mass in the cold gas compo nent also matches w ell the ob¬ 
servational data presented in lPeeples et al.l (l2014l i. We note 
that increasing az can increase the metal mass in the hot 
gas to make the model prediction agree with the data bet¬ 
ter, but it will lower the metal mass in the cold gas com¬ 
ponent so that the agreement with the observational results 
does not hold anymore. Overall, this model nicely repro¬ 
duces the metal masses as functions of galaxy stellar mass 
in all components, except the ejected component, for which 
data constraint is not available. 

Model-RI predicts a trend of the metal mass in stars 
and cold gas that is similar to that in Model-EJ due to the 
common assumption of ejective feedback. However, the two 
models predict very different trends for the metal mass in the 
hot gas and the ejected gas. Since the reincorp oration is fast 
for high-mass galaxies (iHenriaues et al.ll2013l l. a large frac¬ 
tion of metals can return back to the hot halo. Consequently, 
the predicted metal mass in the ejected component first in¬ 
creases with increasing stellar mass, and then decreases, for 
galaxies with stellar masses > 10^° Mq. The metal mass in 
the halo gas predicted by this model increases rapidly with 
increasing galaxy mass owing to both the increase of total 
metal production and the halo mass-dependent reincorpo¬ 
ration. Furthermore, since a large fraction of metals stays 
in the halo, the model predicts high metal mass in the cold 
phase, especially for low-mass halos where cooling from the 
halo gas is efficient. 

Model-PR adopts a moderate ejection with a mass¬ 
loading factor equal to unity, and Fig. [5] shows that even 
with such a small loading factor the ejected mass can still 
contain a large fraction of all the produced metals. Most of 
the metals that are not ejected are predicted to be locked 
in stars and cold disk gas. As in Model-RI, metals in the 
hot halo gas are due to the assumed direct deposition in the 
halo. Again, we see that for the particular choice of az = 0.1, 
the model predicts about 10^ Mq of metals in the halo of 
a galaxy with M* = 10^'^ Mq, which is consistent with the 
observational results. The predicted metal mass in cold gas 
matches the observation result, but the model seems to over¬ 
predict the metal mass in stars. As we will discuss in the 
following section, this discrepancy is the only main problem 
for this model, and suggests that the assumption of uniform 
mixing of metals in the ISM may be invalid. 

The relation between the gas-phase metallicity and 
galaxy stellar mass has been studied observationally 
over the redshi ft rang e 0 < z < 3 b y many authors 
( [Tremonti et al.l 120041 : lErb et ^ l2006l : iMaiolino et al.l 
I2OO8I . e.g.f. However, hydrodynamic simulations and semi- 
analytic models of galaxy formation have difficulty in 


reproducing the relation over the redshift range 

orobed 

bv the observation (iDe Lucia. Kauffmann & White 

l2004 

de Rossi. Tissera & Scannaoiecd 2007: Mouhcine et al.l 

2008|). In Figure [ 6 l we show the predictions of our 

three 


models for the oxygen metallicity in cold disk gas at 
three different redshifts, z = 0, 1, and 2. For compari¬ 
son, we plot three sets of observational results obtained 


by Andrews fc Martinil (l2013lj. iMaioIino et alJ ll2008ll. 
IZabM et al.l ( 2013l f . Yabe et ahl ( 2012l l and lErb et ^ 
( 2006tl . The ejective feedback model (Model-EJ) predicts 
only a weak evolution in the oxygen metallicity - stellar 
mass relation from z = 2 to the present day, which is 
at odd with the increase of the metallicity with time 
for a fixed stellar mass shown in the observational data 
([Mjmnuccj_^tQai][2009UMa et al.l[^08l:lMannucci et al.l 

I 2 OIOI : IZahid et al.ll2013l l. The predicted metallicities at low 
redshifts are too low to match the observation. The reason 
for this is twofold. First, the model over-predicts the cold 
gas fraction for a given stellar mass, which leads to a 
reduced metallicity for the gas. Second, metals are also 
ejected out of galaxies with outflow, further reducing the 
metallicity. The under-predi ction and the la c k of ev olution 
also appear in Model-IR. IHenriaues et al.l (l2013l f found 
that their reincorporation model can match the stellar 
metallicity - stellar mass relation for local galaxies, as¬ 
suming a metal yield p = 0.047 that is 60% higher than 
what is assumed here. However, our results show that the 
reincorporation model is not able to reproduce the evolution 
in the gas phase metallicity unless p is time-dependent. The 
preheating model Model-PR predicts a moderate level of 
evolution and a large scatter in the gas phase metallicity for 
given stellar mass, and the model predictions are consistent 
with the observational data assuming p = 0.03. 

Since we assume that metals in the halo gas is perfectly 
mixed in a density profile, we can make predictions for the 
column density of any chemical species as a function of the 
impact parameter (distance from galaxy center). The COS- 
Halos survey has measured the OVI column density profile 
for a sample of L* galaxies around redshift z = 0.4. Here 
we make the corresponding predictions for a 10^^ Mq halo 
at z = 0.4 and compare our mo del predictions w i th the 
data. In the calculation, we follow iTumlinson et al] (1201 il l 
and assume that the ionization fraction of OVI, fovi = 0.2. 
As shown in Fig[3 all the three models can roughly repro¬ 
duce the shape of the column density profile. The predic¬ 
tions of Model-EJ and Model-PR agree with the data rea¬ 
sonably well, but the column density profile predicted by 
Model-RI is about 2 orders of magnitude too high. This 
over-prediction by Model-RI is again due to the reincorpo¬ 
ration of the ejected gas, which is enriched in metal. 

IWerk et al.l (l2014l l have estimated the mass and column 
density profile of cool gas in the circum-galactic medium for 
44 Milky Way sized galaxies in the COS-Halos sample. They 
found that the total amount of cool gas (with temperature 
in the range 10^ K < T < 10® K) is larger than about 6.5 x 
10^° Mq , and that the cool gas column density profile can be 
well described by a shallow power-law with a power index 
—1 or lower. Assuming the cooling halo gas in our model 
to be this cool component in the circum-galactic medium, 
we compare our model predictions for a typical 10^^ Mq 
halo at 2 = 0.2 with the observational results. To do this, 
we estimate the amount of the halo gas at different radii 
that can cool in a given time interval based on the profile 
of the hot halo gas. The gas cools from a shell at radius r; 
covers a 47rrf area. Assuming that the cooled gas falls onto 
the halo center in one free-fall time from the radius with a 
constant speed, the cool halo gas that drops from radius rt 
has a radial distribution proportional to where r < ri. 
Thus, the total mass density of the cooled halo gas from all 
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Figure 5. The relation between metal mass in different baryonic mass components and the stellar mass for model galaxies predicted by 
different models. The model name is noted in each panel. The red dots denote metal mass locked in stars, the blue dots denote metal 
mass in cold gas in galaxies, and orange dots denote metal mass in hot halo, and the purple dots denote metal mass in gas that is ejected 
out of the halos. The horizontal orange bar marks the observational estimate of the metal mass in halo gas around galaxi es with masses 
over the range of stellar masses covered by the bar as measured in COS-Halos through OVI line dTumlinson et al]|201lh . Solar oxygen 
abundance pattern is assumed. 
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Figure 6. The metallicity of the cold gas (ISM) as a function of stellar mass at z = 0, 0.5, 1, a nd 2 predicted by Model-E J, Model-IR, and 
Mode l- PR. The fitting fun c tions of the mean relati ons derived from the observational data of lAndrews fc Martinil ll2013lh FMaiolino et alj 
ll2008lj . [Zahid et al.l ll2013l‘ ). [Vabe et al.l ll2012f) and lErb et al.l ll2006l^ are shown in each corresponding panel for comparison. 
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Figure 7. The column density of OVI as a function of impact parameter (projected radius from galaxy center) for a 10^^ Mq halo 
predicted by different models. The blue data points are observational data for star forming galaxies and the red data points are for 
quiescent galaxies obtained in iTumlinson et al.fll2011^ ■ 
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Figure 8. The cooling gas column density as a function of impact parameter (projected radius from galaxy center) for a 10^^ Mq halo 
predicted by diffe rent models. The green and orange bands represent the ‘preferred’ and lower limits of the observational estimates of 
IWerk et al.l ll2014|j . 


the outer shells at a radius r can be written as 

Pcooi(r) = ^ y A(ri)rft(ri)ridri, (8) 


where A(ri) is the cooling rate of the halo gas per unit vol¬ 
ume at Ti , and rg is the free-fall time of the gas element from 
the radius where it cools. The free-fall time is estimated as 


Tff = 14.4Gyr 



'M{< n)' 

V IMpc J 

1012 M© 


(9) 


where M(< Vi) is the total mass (dark matter and baryon) 
within radius ri. We then calculate the column density from 
the radial profile and the cool gas column density profiles 
predicted by the three different models, which are shown in 
Fig. [8] The amplitudes of the column density profiles pre¬ 
dicted by both Model-EJ and Model-PR agree with the ob¬ 
servational results, while that predicted by Model-IR seems 


too high. These are consistent with the predicted cooling 
rates shown in Fig.[2] Model-EJ and Model-PR predict a 
much lower rate than Model-IR. More interestingly, both 
ejective feedback models predict profiles that are signifi¬ 
cantly steeper than the observational estimates. In contrast, 
the preheating model predicts an extended cooling gas pro¬ 
file similar to the observational results. 


5 DISCUSSION AND SUMMARY 

Using a set of semi-analytic models of galaxy formation, 
we have demonstrated that the gas and metal contents of 
the galaxy population provide important constraints on the 
feedback processes assumed in galaxy formation models. 
In particular, the masses and the metallicities of baryonic 
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Mass Scale 

Quantity 

Model-EJ 

Model-RI 

Model-PR 

Figure 

Mvir = 10 ^^ Mq 

* 

II 

o 

OK 

too low 

OK 

1 


Mcoid(z=0) 

too high 

too high 

OK 

1 


SFH 

too low at low-z 

too low at low-z 

OK 

2 


/cold(^ = 0) 

too high 

too high 

OK 

3 

Mvu = 1012 M0 

M,{z = 0 ) 

OK 

OK 

OK 

1 


Mcoid(z=0) 

OK 

OK 

OK 

1 


SFH 

OK 

OK 

OK 

2 


/cold(^ = 0) 

OK 

OK 

OK 

3 

M* = 10 ® - 10 i° ® Mq 

Zism(^ = 0) 

too low 

too low 

OK 

6 


^ism(^ = 2 ) 

OK 

OK 

OK 

6 


Mz,,(2 = 0) 

OK 

OK 

too high 

5 


Mz.hot 

OK 

too high 

OK 

5,7 


Table 1. A summary of the success and failure of the three models based on comparisons 
with different observational data. The first column lists the masses of model halos or 
galaxies. The second column lists the quantities in the comparison. The third to fifth 
columns list the performances of the three models. The last column lists the corresponding 
figures from which the conclusion is drawn. 


matter in different phases can be used to distinguish two 
broad classes of feedback scenario, namely ejective and pre¬ 
ventative feedback. We have found that when the models 
are tuned to predict a correct stellar mass for a given halo 
mass, different feedback mechanisms work differently in re¬ 
distributing the metal-enriched material, and hence predict 
very different trends for metallicities. 

We have considered three representative models. Two 
of them implement ejective feedback, in which halos accrete 
baryonic matter together with dark matter at the cosmic 
baryon fraction. To match the observed stellar mass fraction 
of galaxies in the ejective scenario, a significant fraction of 
baryonic matter that cools into the central galaxy has to be 
ejected out of the galaxy by star formation feedback. One 
of the ejective models implements a simple assumption that 
the ejected material never comes back, while the other one 
considers a more sophisticated model in which the ejected 
gas can be re-incorporated into the halo. The third model 
implements a preventative feedback model, in which feed¬ 
back works to preheat the IGM, and the thermal pressure 
prevents a fraction of baryon from collapsing into low-mass 
halos. To explore the characteristics of these models, we have 
contrasted their predictions in the masses and metallicities 
of stars, cold disk gas, hot halo gas, and ejected gas. Table 
1 summarizes the success and failure of the model predic¬ 
tions in comparison with existing observational data. In the 
following we discuss the implications of these results. 

First, ejective and preventative models predict different 
fractions of baryons that reside in a halo and outside the halo 
viral radius. The pure ejective model predicts that most of 
the baryon mass is ejected out of the halo, and only a small 
fraction is retained in the central galaxy and the halo. The 
ejective model with gas reincorporation has similar predic¬ 
tions for low-mass halos, but for Milky Way sized galaxies 
it predicts that the dominating fraction of baryons is in the 
halo because the ejected gas can return quickly back to the 
halo. 

Second, we have shown that when the outflow mass¬ 
loading parameter is tuned to be large enough in the ejective 
models so that the predicted total stellar mass matches the 


observed stellar mass - halo mass relation, model galaxies 
tend to have too high a gas mass fraction compared with 
observational data. Meanwhile, the ejective models also pre¬ 
dict a too low gas-phase metallicity. Although stronger out¬ 
flow with an even larger mass-loading factor could reduce 
the gas fraction in the ejective models, the metal mass will 
inevitably decrease further, resulting in an even lower gas- 
phase metallicity. Thus, the ejective models appear to be 
unable to accommodate the observed gas/stellar mass frac¬ 
tion and gas-phase metallicity simultaneously. 

Third, regardless of the details of how ejection and 
reincorporation are parameterized, the ejective models al¬ 
ways predict a gas phase metallicity-stellar mass relation 
with roughly fixed amplitude and slope since z = 2 without 
strong time evolution, which is inconsistent with current ob¬ 
servations. 

Finally, we have shown that if feedback is primarily 
ejective and the ejected mass can be reincorporated into 
massive halos rapidly as suggested in recent models (e.g. 
iHenriaues et al.|l2013h . gaseous halos may contain too much 
metal mass to be consistent with COS-Halos observational 
results. 

In contrast the preventative model can accommodate 
almost all observational data, including the evolution of 
masses and metallicities of baryonic matter in different 
phases. The only problem is that the model predicts too high 
metal masses in the stelar component of low-mass galaxies. 
Given that this model predicts the correct metal mass in 
ISM, and the correct total stellar and cold gas masses, the 
problem cannot be solved by assuming stronger outflows 
with uniform mixing. The over-prediction suggests that a 
significant fraction of stars in low-mass galaxies may form 
in a medium with a significant lower metallicity than the 
average ISM metallicity. This can happen if a fraction of 
star formation in low-mass galaxies occurs in early starburst 
before z = 2, when the ISM metall icity has not been sub¬ 
stantially enriched iILu et al.ll2014bh . Alternatively, a large 
fraction (> 70%) of the stellar mass in low-mass galaxies 
may have formed in regions where the metallicity is lower 
than the average. Indeed, as demonstrated in detail in Lu et 
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al. (2015, submitted), both the observed stellar-phase metal- 
licity and gas-phase metallicity as functions of galaxy stellar 
mass can be reproduced, if star formation occurs in an ISM 
with inhomogeneous metallicity, where a large fraction of 
the star forming ISM has metallicities below the average 
ISM metallicity. Since real ISM i s inhomogeneou s in metal¬ 
licity (ISmartt fc RollestonI IiOOtI : IHo et alTl2015h . this only 
discrepancy for the preventative model may not be a real 
problem after all. 

In summary, the results we have obtained indicate that 
the basic assumption of the ejective model that all bary- 
onic matter first collapses into a halo and the baryonic mat¬ 
ter budget in a galaxy is regulated by outflows is problem¬ 
atic. This, together wit h the success of the preven tive model 
considered here and in ILu. Mo fc Wechsleif ll2015l ) , suggests 
that feedback is likely preventative in nature. In the present 
paper we have explored the different models by manually 
varying the key model parameters. These parameters cer¬ 
tainly have degeneracy between them. In addition, the two 
processes, ejection and prevention, are expected to be de¬ 
generate to some extent for some of the observational con¬ 
straints. It remains to be revealed how these two processes 
are degenerate and what observational data can effectively 
break the degeneracy. We will investigate these questions 
using the Baye sian model inference approach developed by 
ILu et all (l201ll ') in a future paper. 
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